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It is shown that at low temperatures and moderate electron dephasing the electron transmission 
function reveales a structure containing information about donor/acceptor sites effectively partic- 
ipating in the electron transfer processes and primary pathways of electrons tunneling through 
molecular bridges in macromolecules. This important information can be obtained as a result of 
analysis of experimental low temperature current-voltage characteristics for chosen molecules. 
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It has been estabilished that molecular electron trans- 
fer (ET) is essentially a combination of nuclear environ- 
ment fluctuations and electron tunneling. Due to the 
large distances between donor and acceptor, ET is mostly 
provided by intervening molecular bridges, giving rise to 
a set of intermediate states for the electron tunneling . 
The expression for the ET rate including both electronic 
and nuclear factors was first proposed by Marcus j^l^Q 
and can be written as follows: 

KeT = KelKnVn ■ (1) 

Here, Kf,i is the electron transmission coefficient, Kn is 
the nuclear transition or Franck-Condon factor, and Vn 
is the effective nuclear vibration frequency. 

The enormous size of biological molecules makes calcu- 
lations of if extremely difficult In 
this paper a simple approach is proposed which enables 
us to show that some intrinsic characteristics of the in- 
tramolecular ET such as pathways of tunneling electrons 
and donor/acceptor coupling to the bridge at different 
values of tunnel energy E can be obtained in experiments 
on low-temperature electrical transport through the cor- 
responding molecules. These data niay be available under 
moderate electronic dephasing [ij, jlj, [iJI . A similar ap- 
proach could also be developed to analyze characteristics 
of intermolecular ET processes in extended biosystems. 

Existing theories of electronic transport through 
macromolecules are mostly based on a very simple model 
simulating both electrodes by semi-infinite tight-binding 
chains attached to the ends of the molecule which is 
also simulated as a single tight-binding chain of sites 
[H m Q E H [13. Here we adopt a resembling 
model. However, we take into account that donor and 
acceptor subsystems in realistic macromolecules are com- 
plex and include many sites providing effective coupling 
to the bridge, as well as to the attached electrodes. 

To simplify further calculations we assume coherent 
electron tunneling through the macromolecule to be the 
predominant transport mechanism. Correspondingly we 
treat the electron transport through the molecule as a 
sequence of tunnelings between potential wells. Each 



well represents one atomic orbital of the donor-bridge- 
acceptor system. Within this approach, any atom is rep- 
resented by a set of sites corresponding to its states. For 
further simplification we assume than intraatomic hop- 
ping integrals are smaller that interatomic ones, therefore 
we consider electron transfer between different sitesas its 
tunneling between different atoms. 

Assume that sites " i" are associated with those atoms 
of the bridge which have a strong coupling to the donor, 
and sites "j" are associated with atoms interacting with 
the acceptor. Then we can treat any site "i" as the ori- 
gin of a chain along which an electron moves from donor 
to acceptor, and any site "j" as the end of the chain. 
Thus, we have a set of chains (pathways) for an elec- 
tron moving along the bridge. In further analysis, we 
concentrate on the case when the chains weakly inter- 
act so we can consider them separately. We simulate the 
donor /acceptor subsystems as sets of semi-infinite tight- 
binding homogeneous chains. Each chain is attached to a 
site of donor/acceptor which can be effectively coupled to 
the bridge. As before we assume that the chains do not 
interact with each other. The adopted model does not 
enable us to carry out a proper treatment of coupling 
of electrodes to the molecule at metal-molecule-metal 
junctions which is a nontrivial problem [lil Hol Eol EH. 
However, it seems reasonable to conjecture that in the 
molecules with complex donor and acceptor subsystems 
this process do not strongly depend on details of coupling 
of the electrodes to the donor/acceptor due to a com- 
paratively large size and complicated structure of these 
subsystems. 

In the following calculations we start from an effective 
tight-binding Hamiltonian for a single chain included into 
the bridge: 



H,jj = Ha + Hi+Hv + HA. (2) 

Here first two terms describe the chain itself, and their 
matrix elements between states \k > and \l > correspond- 
ing to the fcth and Ith sites are given by: 
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iHo)ki = ak6, 



kl 



{Hi)ki — Vki 



(3) 



where Vki= when k = I and only states associated with 
valence electrons are considered. The diagonal matrix 
elements ak are ionization energies of electrons at sites 
k, while Vki — Vik includes both direct and exchange en- 
ergy contributions for an electron to transfer between the 
kth and ^th sites. Remainig terms represent self-energy 
corrections arising due to the coupling of the donor (-ffp) 
and acceptor (i/^) to the bridge: 
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E ~ tn — (Jn 



(4) 
(5) 



Here, 2?mi and Ajn are, respectively, coupling strengths 
between the mth donor site or the nth acceptor site 
and the «th or jth site of the bridge, and cr„j_„ = 

\ \Ora,n " i^J'^lm^n " ^'m,n} are the sclf-cnergy correc- 
tions of the semi-infinite chains attached to the corre- 
sponding sites ,12]. The parameters ^m^„ — E — tm.m 
where Cm.n and 7m^„ are the ionization energies of elec- 
trons at the corresponding donor/acceptor sites, and 
the nearest- neighbor hopping integrals for the chains. 
Summation in Eqs. (4) and (5) is carried out over all 
donor/acceptor sites coupled to the bridge. Due to the 
presence of the self-energy corrections the eigenvalues of 
the effective Hamiltonian [Eq. (2)] include imaginary 
parts which give broadening of the bridge energy lev- 
els Ei. The energy levels are broadened further if include 
scattering processes in the bridge which are not consid- 
ered here. 

An electric tunneling current / flowing from donor to 
acceptor through the bridge in the presence of a small 
voltage V applied across the bridge has the form |^ 

e f°° 

/=— / dE T{E)[f {E - ^Ji^) - f {E - , (6) 

where f{E) is the Fermi function, and the chemical po- 
tentials fii and fi2 are determined by the equilibrium 
Fermi energy of the bridge Ep and the effective voltage 
V across the bridge 0, i.e. 

^1 = i?F + (1 - v) eV"; fi2 = Ep - rjeV. 

Here, the parameter rj characterizes how the voltage V is 
divided between the two ends of the bridge. The electron 
transmission function is given by the formula 



nE)^2Y,^^\Gvf^J 



(7) 



The summation in Eq. (7) is carried over the bridge 
states \i > and \j > therefore contributions from all pos- 
sible pathways are contained here. The quantities A^j 
are imaginary parts of the self-energy corrections 
and Gij is the matrix element of the Green's function 
corresponding to He/ / : 



{^\iE 



Ho — Hi — Hj) — Hj) 



(8) 



We see that the dependence of the electron transmission 
function on energy is determined by the contributions 
from different donor/acceptor sites, as well as from the 
Green's function matrix elements corresponding to dif- 
ferent chains included into the bridge subsystem. 

We can easily arrive at the first approximation for Gij 
neglecting the broadening of the bridge energy levels. 
Taking into account only nearest neighbors (NN) and 
next nearest neighbor (NNN) couplings, we can intro- 
duce a new notation for the nonzero hopping integrals, 
i.e., Vk,k+i = Vk~i,k = Pk and Vk,k+2 = Vk-2,k = Ik- 
Then we can calculate the Green's function matrix ele- 
ments for a given chain which makes a start at the site " i" 
and ends at the site "j", applying a diagram technique 
similar to that described in detail for the (NN) approxi- 
mation 23] . We present our chain of sites as a graph , 
using the following rules: 

(i) Any site included in the chain corresponds to a 
graph vertex. 

(ii) A term E — ak corresponds to a loop attached to 
the vertex "fc" . 

(iii) A hopping integral Vki corresponds to an edge orig- 
inating at the vertex "fc" and ending at the vertex . 

As a result we obtain the graph shown below: 



ak E - ak+i E 




Ij 7j + l 7j+2 7fc-2 7/c-l Ik Ik+l lk+2 

E — oij-2 E — aj^i E — aj 




lj-4 7j-3 lj-2 



Now we consider all possible cycles for this graph and 
define the value of a cycle O as the product of values 
of the edges included into the cycle multiplied by the 
number of the edges connecting the firstand last sites 
with a negative sign. When the numbers of the edges 
associated with the paths from the first site to the last one 
and back differ, we multiply by the maximum number of 
edges. When the cycle includes only one vertex, its value 
is defined as E — ak where k is the number of the vertex. 
Then we introduce a cyclic term which is a set of cycles 
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including all sites of the system. The value of the cyclic 
term is equal to the product of the cyclic values of these 
cycles. The sum of all possible cyclic terms O,- gives us 
the cyclic value 6(r,y) of the graph Fy. An expression 
for the matrix element Ga can be obtained in the form: 



(9) 



Here, the summation is carried out over all possible path- 
ways along the chosen chain. The value of P/J equals the 
product of all edges along the pathway. When the path- 
way includes all vertices of the graph we have 



(10) 



k—i 



When the pathway does not include some vertices, this 
produces different values for the corresponding . For 
example, when the electron isbeing transferred from the 
site I to the site 1 + 2, passing the intermediate site Z -I- 1, 
the corresponding edge is equal to {E — ai+i)^i. Conse- 
quently, we obtain for the pathway which excludes this 
vertex: 



P 

J J.I 



{E - Q;i+i)7/ 



n 



f3k 



(11) 



Following a similar way we can get explicit expressions 
for all terms in the numerator of Eq. (9). The expres- 
sion (9) can be significantly simplified if we assume that 
(NNN) couplings are small compared to (NN) couplings. 
Then we obtain the following approximation: 



1 



i-i i-2 j-i 

l=i k=i;k^Ll+l\ 



\k=l 



(12) 

where 9o(Fy ) is the cyclic value of the graph calculated 
within the (NN) approximation. 

The Green's function matrix elements [Eq. (9)] have 
poles which can be detemined by solving the equations 



e.,(F) = o. 



(13) 



These poles correspond to energy states of a given chain. 
The set of these energy states characterizes the chosen 
chain (pathway) and different chains have different sets of 
energy states. When we take into account the broadening 
of the bridge energy levels, the singularities of the Green's 
function matrix elements are converted to the peaks of 
finite height and width. 

To apply the described method of calculation of the 
matrix elements for a real macromolecule we have 
to determine which atomic orbitals have to be used to 
arrange a chain between the sites "i" and "j". To choose 
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FIG. 1: The calculated current (nA) - voltage (V) charac- 
teristics for different temperatures. 



suitable orbitals we can employ a rigorous quantum me- 
chanical method of tunneling interatomic currents [2^ 
which was alredy successfully used to separate out elec- 
tron transfer pathways in long-range ET processes [isf . 
Data concerning the electronic structure of a considered 
molecule can be obtained within an extended Huckel 
or some other approximation 26] . As a result we get the 
set of chains which represent the bridge, and we can use 
our result [Eq. (9)] in further analysis. 

It follows from the obtained results [Eqs. (7) and 
(9)] that the electron transmission function contains 
important information concerning intramolecular elec- 
tron transfer. Under low temperatures and moderate 
electronic dephasing T{E) can exhibit scries of peaks. 
Their location is determined with contributions from 
donor/acceptor sites predominating in the ET at a given 
interval of tunneling energy E, and with the energy spec- 
trum of a chain of sites connecting them. When E takes 
on a value close to ek,i or to one of the poles of the Green's 
function the corresponding term in Eq. (7) can surpass 
all remaining terms. Thus, for different values of E dif- 
ferent donor/acceptor sites and different pathways can 
predominate in the electron transfer through the bridge. 

The structure of T{E) can be revealed in the low 
temperature current-voltage characteristics for the elec- 
tronic transport through the molecular bridges [Eq. (6)]. 
As a result of numerical simulations of electrical trans- 
port throu gh p orphyrin-nitrobenzene macromolecules it 
was shown |27l| that at low temperatures current-voltage 
curves can exhibit a steplike behavior which disappears 
as the temperature raises (see Fig. 1). This originates 
from a steplike character of Fermi distribution functions 
at low temperatures along with the "comblike" structure 
of T{E). At a given voltage V the difference of Fermi 
functions in Eq. (6) takes on nonzero values only in the 
interior of a certain energy range including Ep. There- 
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fore the magnitude of the current / at a given voltage 
is determined by the contributions of peaks of T{E) lo- 
cated in this energy range. When the applied voltage 
increases, this enhances the width of the relevant energy 
interval. The tunneling current / abruptly changes when 
an extra peak of the electron transmission function ap- 
pears there. Widths of the plateaus are equal to the 
distances between adjacent peaks of T{E) and magni- 
tudes of sudden changes in the current correspond to the 
heights and shape of these peaks. At higher temperatures 
Fermi functions lose their steplike character and the pla- 
teous are washed out. Another reason for the structure 
of T{E) (and I — V curves) to be eroded is the electronic 
phase-breaking effect which arises due to stochastic fluc- 
tuations of the ion potential. We conjecture, however, 
that at low temperatures the electronic dephasing effects 
could be reduced so that the structure of T{E) could be 
revealed. 

On these grounds we believe that experiments on elec- 
tronic transport through macromolecules at low tem- 
peratures (T ~ IK) could provide an additional in- 
formation upon characteristics of the electron transfer. 
Namely, comparing the structure of the electron trans- 
mission function reconstructed on the basis of experi- 
mental I — V curves with that obtained as a result of 
calculations using Eqs. (7) and (9), we can make conclu- 
sions concerning actual primary pathways of electrons, as 
well as sites of the donor/acceptor subsystems involved 
in the ET process at different values of the tunneling en- 
ergy. The model adopted in the present work provides 
us with the results suitable for a quantitative compari- 
son with the results of proposed experiments even avoid- 
ing a proper and reliable calculation of the equilibrium 
Fermi energy of the bridge, and the effective voltage V . 
The lack of information about proper values of Ep and 
V produces an uncertainty in the location of the origin 
at the "voltage" axis in the Fig. 1, so we cannot iden- 
tify steps of the I — V curves separately. Nevertheless, 
changes in the values oi Ep and V do not influence the 
electron transmission function, therefore the structure of 
series of the peaks remains fixed. This enables us to iden- 
tify some series of peaks analyzing sequences of widths 
of successive steps of the current- voltage curves. Such 
analysis also can give reasonable estimations for the Ep 
and V for the chosen molecule [2^ . 
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